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ABSTRACT 

About 10% of the massive main sequence stars have recently been found to host a strong, large scale magnetic field. Both, the origin 
and the evolutionary consequences of these fields are largely unknown.We argue that these fields may be sufficiently strong in the 
deep interior of the stars to suppress convection near the outer edge of their convective core. We performed parametrised stellar 
evolution calculations and assumed a reduced size of the convective core for stars in the mass range 16Mo to 28 Mq from the zero 
age main sequence until core carhon depletion.We find that such models avoid the coolest part of the main sequence hand, which 
is usually filled by evolutionary models that include convective core overshooting. Furthermore, our ‘magnetic’ models populate the 
blue supergiant region during core helium burning, i.e., the post-main sequence gap left by ordinary single star models, and some of 
them end their life in a position near that of the progenitor of Supernova 1987A in the HR diagram. Further effects include a strongly 
reduced luminosity during the red supergiant stage, and downward shift of the limiting initial mass for white dwarf and neutron star 
formation. 
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1. Introduction 


Massive stars are important cosmic engines, because they drive 
the evolution of star forming galaxies (Mac Low et al. 2005), 
and are very bright and thus visible from very long distances 
(Kudritzki et al. 2008). Nevertheless, even their longest lasting 
evolutionary stage, core hydrogen burning, has not yet been un¬ 
derstood well. While one of the two culprits, mass loss, strongly 
affects only the very massive stars, it is the other one, internal 
mixing processes, that still leads to many questions. This refers 
to thermally driven mixing as caused by convection, overshoot¬ 
ing, or semiconvection as well as to rotationally induced mixing 
and the mixing affected by magnetic fields. 

In recent years, it has become evident that the concept of 
considering one evolutionary path for stars of a given mass fails 
in the massive star regime. Hunter et al. (2008) and Brott et al. 
(2011) show that early B-type main sequence stars fall into sev¬ 
eral distinctly different categories, as classified by their spin and 
surface nitrogen abundance. Dufton et al. (2013) show that the 
rotational velocity distribution of EMC early B stars is bimodal, 
similar to late B and A stars (Royer et al. 2007). The physical 
reasons for this are not clear yet, but could relate to initial rota¬ 
tion, binarity, and strong internal magnetic fields (Langer 2012). 

In the present paper, we consider one of the least investi¬ 
gated aspects of massive star evolution, namely the effects of 
internal magnetic fields. Four different types of magnetic fields 
in massive stars have been discussed in recent years. First, there 
is the idea that differential rotation in their radiative envelopes 
produces toroidal fields through a dynamo process (Spruit 2002). 
While the functioning of this dynamo is still being debated (Zahn 

2011) , asteroseismic measurements in red giants (Mosser et al. 

2012) call for a strong angular momentum transport mechanism 
which could be provided by this predominantly toroidal field 
(Suijs et al. 2008; Heger et al. 2005). This type of field, which 
would not be directly observable, and its consequences would be 
as ubiquitous as differential rotation. 


Second, massive stars have convective cores, and it appears 
likely that an afi-dynamo can produce a magnetic field in these 
cores. In fact, the MHD simulations of Brun et al. (2005) show 
the functioning of such a dynamo in the convective cores of even 
slowly rotating 2 M© stars, where the field growths saturate near 
their equipartition values. MacGregor & Cassinelli (2003) show 
that this field is unlikely to rise to the surface of the star during 
its lifetime, so again, it may not be directly observable. Its effects 
in the vicinity of the convective core are largely unexplored. 

A third type of magnetic field in massive stars has been 
suggested by Cantiello et al. (2009), who propose that an aQ- 
dynamo may work in their iron opacity driven sub-surface con¬ 
vection zones. The produced flux tubes could then buoyantly 
float to the stellar surface and produce magnetic spots (Cantiello 
& Braithwaite 2011), which may be observationally indicated by 
the moving discrete absorption components in the UV resonance 
lines found in O-star spectra (Prinja & Howarth 1988; Kaper et 
al. 1997). 

Here we are concerned with a fourth type of magnetic field in 
massive stars, namely strong large scale fields that are anchored 
in the radiative stellar envelope and that form an observable mag¬ 
netosphere outside the star. Only OB stars with such fields will 
be designated as magnetic OB stars in this paper. Until recently, 
only three magnetic O-type stars were known, HD 37022 (Do- 
nati et al. 2002), HD 191612 (Donati et al. 2006) and f Orio- 
nis A (Bouret et al. 2008). However, new intense searches have 
revealed many more magnetic OB stars (Grunhut et al. 2009, 
Petit et al. 2013, Hubrig et al. 2013, 2014), thus implying an 
incidence fraction of strong, large scale B-fields in OB stars of 
about ten percent (Grunhut & Wade 2012). 

In contrast to the first three types of magnetic fields in OB 
stars, which all need a dynamo process to continuously produce 
the field because it would quickly decay otherwise, the strong 
large scale fields are thought to be stable for two reasons. First, 
they occur mostly in very slow rotators and in stars with radia- 
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tive envelopes, such that there is no known dynamo process that 
could produce them. Second, Braithwate and Spruit (2004) find 
that magnetic field configurations resembling the ones observed 
in magnetic OB stars can indeed be stable, in the sense that their 
decay may take longer than the stellar lifetime. 

The origin of these stable, large-scale magnetic fields in mas¬ 
sive stars is mainly described by two ideas. The hypothesis of 
fossil fields (Braithwaite & Nordlund 2006) assumes that mag¬ 
netic fields in a molecular cloud are conserved during star forma¬ 
tion. It is, however, difficult to explain, why only a low percent¬ 
age of stars preserves a magnetic field in this formation channel. 
The effects of a fossil magnetic field on convective core dynamos 
in A-type stars are described in Featherstone et al. (2009). 

Magnetic stars can also be the result of strong binary inter¬ 
action (Langer 2014). For example, they could be the product 
of the merger of two main sequence stars, which occurs quite 
frequently because a large number of massive stars have close 
companion stars (Sana et al. 2012). Garcia-Berro et al. (2012) 
show that strong magnetic fields can form in the hot, convec¬ 
tive, differentially rotating envelope of a merger product of two 
white dwarfs, which may occur similarly in the case of a main 
sequence merger. Calculations by de Mink et al. (2014) suggest 
that among the galactic early-type stars, 8)^^ % are the products 
of a merger event from a close binary system. The order of mag¬ 
nitude obtained in this study agrees well with the occurrence of 
magnetic massive stars. Furthermore, that many magnetic OB 
stars are nitrogen-enriched at their surface complies well with 
the models of merger stars by Glebbeek et al. (2013). 

The evolutionary consequences of strong large scale mag¬ 
netic fields in the stellar interior are largely unexplored. Most 
observed magnetic stars rotate slowly, such that rotational mix¬ 
ing plays no role at least at the present time. Furthermore, it is 
found that in many situations, large scale magnetic fields and 
convection repulse each other (e.g., Gough & Tayler 1966; Gal¬ 
loway & Weiss 1981). Therefore, while in evolutionary models 
of massive stars it is usually assumed that the convective core 
is extended over its standard size during core hydrogen burning 
(Brott et al. 2011; Ekstrom et al. 2012), we assume the oppo¬ 
site in the present paper, which means that internal large scale 
magnetic fields restrict the convective core to a smaller region 
than what the Ledoux criterion for convection would predict. 
We elaborate our assumptions and the computational method in 
Sect. 2, and present our evolutionary calculations in Sect. 3. Af¬ 
ter a comparison with observations in Sect. 4 we give our con¬ 
clusions in Sect. 5. 

2. Method 

Our stellar models are calculated with a hydrodynamic stellar 
evolution code, which is described in detail in Yoon et al. (2006) 
and Kohler et al. (2015). 

Convection is modelled using the Ledoux criterion, assuming a 
mixing-length parameter 1-5 (Langer 1991). Semiconvec¬ 
tion is implemented as described by Langer et al. (1983), apply¬ 
ing a value of asem = 0.01 (Langer 1991; see also Zaussinger 
and Spruit 2013). The description for mass loss follows Vink 
et al. (2000) for main sequence stars and blue supergiants, and 
Nieuwenhuijzen et al. (1990) for red supergiants. Lor our models 
we adopt solar abundances as the initial composition as in Brott 
etal. (2011). 

The observed fields in magnetic massive stars are large scale 
fields, and while some of them have a more complex geome¬ 
try many have a dominant dipole component (Donati & Land- 
street 2009). Braithwaite & Spruit (2004) show that in order 


to achieve a stable magnetic field configuration, a dipole field 
needs to be balanced by an internal toroidal field. While a typi¬ 
cal value for the observable part of the dipole field is 1 kG (Los- 
sati et al. 2015), magnetic flux conservation would then predict 
field strength of up to 10® G near the convective core of a massive 
star. Additionally, Braithwaite (2009) has shown that the toroidal 
component can have a much higher field strength in a stable con¬ 
figuration than the poloidal component. Therefore, field strength 
of the order of 10*^ G cannot be excluded in the deep interior of 
magnetic massive stars and this corresponds to the equipartition 
field strength and indicates direct dynamical effects of the mag¬ 
netic field on the convective motions. 

Additionally, Zahn (2009) has shown that the magnetic field 
strength that is required to suppress convection is about two or¬ 
ders of magnitude less than the equipartition field strength, when 
thermal and Ohmic diffusion are taken into account. Zahn ar¬ 
gues that this effect is observationally confirmed by the pres¬ 
ence of sunspots and by surface inhomogeneities in Ap/Bp stars. 
We therefore assume in the following that convection near the 
boundary of the convective core is suppressed in magnetic mas¬ 
sive stars. 

Lydon & Sofia (1995) have suggested a modification of the 
stellar structure equations to account for the presence of a large 
scale magnetic field inside a star (see also Leiden & Chaboyer 
2012). In our models, we ignore the changes induced by B fields 
on the structure of the star, which implies that we assume that 
the magnetic energy density is lower than the thermal energy 
density. We conduct a simple parametric study of stellar evolu¬ 
tion models here with a diminished core mass compared to what 
is predicted by the Ledoux criterion for convection. To this end, 
we applied a modified Ledoux criterion as 

< f *'^Led, (1) 

where 

^Led = ^ad - —V^. (2) 

XT 

Whereas this mimics the convection criterion to some extent in 
the presence of a vertical magnetic field proposed by Gough & 
Tayler (1966; see also Schatten & Sofia 1981; Tayler 1986), a 
quantitative prediction of the suppressive effects of large scale 
fields upon convection requires multi-dimensional MHD calcu¬ 
lations. We thus refrain from relating our parameter / to mag¬ 
netic field strengths, but make a heuristic choice for its values 
between 1.12 and 1.20, which leads to a reduction of the convec¬ 
tive core mass of the order of 15% in our models. We restricted 
this change to the core hydrogen-burning phase, because the in¬ 
teractions of the magnetic field with the helium-burning core are 
uncertain. In particular, it remains unclear whether a significant 
fraction of the magnetic field will survive below the hydrogen¬ 
burning shell source. 

At present, the origin of the large scale fields in massive stars 
is not clear (Langer 2014). While it might exist throughout core 
hydrogen burning if the field was inherited from the star forma¬ 
tion process, it would be present for a fraction of that time if it 
was created in a binary merger or mass transfer process. How¬ 
ever, a significant mass gain is expected to rejuvenate the star, 
that is to say that the binary interaction resets its evolutionary 
clock. Even though the time is not reset to zero in this case, we 
assume that the effect of a magnetic field decreasing the con¬ 
vective core is unaltered since the beginning of core hydrogen 
burning in our models. 

All our models are evolved through core carbon burning, 
such that their presupernova positions in the HR diagram can 
be reliably predicted. 
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Table 1. Key quantities of the computed evolutionary models. Columns: acronym of the sequence, initial mass Mj, core reduction parameter /, the 
initial mass of the convective core Mcc.i, initial and final helium core masses Mne.i and Mhcj, lifetimes of core hydrogen and core helium burning 
Th and Thc, and presupernova effective temperature and luminosity Lf 


Acronym 

M. 

M© 

/ 

-^cc,i 

Mo 

4^He,i 

M© 

4^He,f 

M© 

Th 

Myr 

THe 

Myr 

Teff.f 

kK 

logiLflLo) 

M16fl.l2 

16 

1.12 

4.8 

2.41 

3.83 

8.75 

1.22 

3.48 

4.65 

M16fl.l6 

16 

1.16 

4.4 

2.27 

3.77 

8.56 

1.36 

3.51 

4.61 

M16fl.l8 

16 

1.18 

4.3 

2.13 

3.61 

8.46 

1.43 

3.54 

4.56 

M16fl.20 

16 

1.20 

4.1 

2.12 

3.52 

8.37 

1.37 

17.4 

4.70 

M18fl.l2 

18 

1.12 

5.7 

2.90 

4.53 

7.56 

0.94 

3.48 

4.73 

M18fl.l6 

18 

1.16 

5.2 

2.68 

4.44 

7.38 

1.04 

3.49 

4.72 

M18fl.l8 

18 

1.18 

5.1 

2.58 

3.73 

7.30 

1.33 

23.2 

4.88 

M18fl.20 

18 

1.20 

5.0 

2.40 

3.32 

7.21 

1.46 

20.5 

4.85 

M20fl.l0 

20 

1.10 

7.0 

3.46 

5.39 

6.78 

0.71 

3.43 

4.88 

M20fl.l3 

20 

1.13 

6.5 

3.28 

4.33 

6.66 

0.92 

22.9 

5.01 

M20fl.l5 

20 

1.15 

6.2 

3.14 

4.15 

6.57 

1.00 

23.4 

5.01 

M22fl.l8 

22 

1.18 

7.0 

3.36 

4.30 

5.82 

0.93 

24.6 

5.10 

M24fl.l8 

24 

1.18 

7.9 

3.86 

5.36 

5.32 

0.65 

25.1 

5.16 

M28fl.l8 

28 

1.18 

10.0 

4.76 

6.54 

4.61 

0.50 

24.7 

5.31 


3. Results 

We performed calculations for stars in the mass range 16 to 28 
Mq, assuming different values for the convective core reduction 
parameter / (cf., Eq. 1), and Table[T| gives an overview of the 
computed sequences. Figure [T] provides Kippenhahn diagrams 
showing the internal evolution of two 16 Mq star as an exam¬ 
ple, where reduction parameters / = 1.12 (top) and / = 1.20 
(bottom) have been adopted. The convective cores are initially 5 
% and 18 % lower in mass than for a non-reduced core, which 
leads to initial helium core masses of 2.4 M© and 2.1 M©, in both 
cases, compared to 2.9 M© for the unreduced case. 

As shown in Fig. [T] extended semiconvective shells develop 
on top of the convective cores for the models with / > 1. While 
the interaction of these semiconvective layers with the large scale 
magnetic field is unclear, it is likely that the field also inhibits 
any semiconvection. For our simulations, this would make little 
difference, since the chosen small semiconvective efficiency pa¬ 
rameter (see Sect. 2) means that there is very little chemical or 
energy transport in these zones. 

The two stars shown in Fig. [1] develop helium cores of ini¬ 
tially about 2.6M© (/ = 1.12) and 2.2M© (/ = 1.20), respec¬ 
tively, which increase to a final mass of about 3.5 M© during 
core helium burning. Such values are expected in stars in the 
mass range IIM© to 12M© for the physics used in the mod¬ 
els of Brott et al. (2011), that is, relatively large convective core 
overshooting. 

Figure |2](top) shows the evolutionary tracks of both 16M© 
models in the HR diagram. The minimum effective temperature 
during core hydrogen burning is higher for the model with the 
larger core reduction parameter, but both are very high com¬ 
pared to that of a model by Brott et al. (2011) shown in the 
same plot. Further on, the tracks in Fig. |2] show that the model 
computed with /=1.20 never evolves into a red supergiant. In¬ 
stead, its lowest effective temperature is about lOOOOK, and it 
ends its evolution at about 17 000K, which is close to the ef¬ 
fective temperature of the progenitor of Supernova 1987A. The 
16M© model computed with /=1.12, on the other hand, moves 
to the red supergiant branch quickly after core hydrogen exhaus¬ 
tion, and remains there for the rest of its evolution. Its luminosity 
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Fig. 1. Kippenhahn diagrams of two 16 Mq models computed with a 
convective core reduction parameter of /= 1.12 (top) and /=1.20 (bot¬ 
tom). Convective (green) and semiconvective (red) regions are shown as 
a function of mass coordinate and time. Thermonuclear energy genera¬ 
tion or loss is colour-coded as shown on the scale on the right-hand side. 
The evolutionary tracks of both models in the HR diagram are shown in 
Fig.[2] 
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during core helium burning and its presupernova luminosity of 
log L/Lq - 4.65 is much lower than the luminosity of models 
that assume convective core overshooting. 

3.1. Variation in the core reduction parameter 

For stars of initially 16 Mq, ISMq, and 20 Mq, we computed a 
series of models with different reduction parameters / (see Table 
[TJ. The evolutionary tracks in the HR diagram of these models 
are shown in Fig.|2] together with tracks for non-rotating Galac¬ 
tic stars from Brott et al. (2011) for 15 Mq and 20 Mq. 

Figure |2]reveals that the tracks in the HR diagram show dis¬ 
continuous behaviour as a function of the parameter /. While 
at 16 Mq, for example, the tracks vary little but in a systematic 
way for / = 1.12, 1.16, and 1.18, with slightly reduced main 
sequence widths and presupernova red supergiant luminosities 
for higher values of /, the track computed with f - 1 .20 be¬ 
haves qualitatively in a different way. Its main sequence part is 
still within the trend of the tracks with lower /-values, but dur¬ 
ing core helium burning, the model avoids the red supergiant 
regime, and it evolves into a presupernova configuration that is 
significantly more luminous than for the three other 16Mo se¬ 
quences. 

Such discontinuous behaviour is also found at ISMq and 
20 Mq, although the jump occurs at different /-values. For the 
18Mo sequences, / = 1.16 still leads to a red supergiant con¬ 
figuration, whereas f - 1.18 and / = 1.20 does not. And in the 
20 Mq models, the discontinuity is found between f — 1.10 and 
/ = 1.13. If / were to measure the strength of the large scale 
magnetic field near the convective core boundary of our models, 
less massive stars would require stronger fields in order to un¬ 
dergo core helium burning as a blue rather than a red supergiant. 

The discontinuous behaviour also concerns the presupernova 
position of the models in the HR diagram. The models that 
burn helium as blue supergiants also end their evolution in this 
regime, while those that burn helium as red supergiants also die 
as such. Whether this remains true for intermediate values of 
/, which are not considered here, is not known. However, the 
16 Mq model computed with / = 0.20 moves rather far towards 
the red supergiant regime after core hydrogen exhaustion and 
then turns around, such that the model burns helium as a red 
supergiant. It might thus be possible that models with /-values 
near the discontinuity have a brief red-supergiant phase before 
turning into core helium-burning red supergiants. Similarly, all 
our core helium-burning blue supergiant models evolve redwards 
during core helium burning. However, those that are still blue 
supergiants at core helium exhaustion will produce blue presu¬ 
pernova stars, since the models evolve mostly bluewards during 
their post-core helium- burning life time. 

Our ‘magnetic’ red supergiant presupernova stars have re¬ 
markably low luminosities, compared to stars computed with 
convective core overshooting, and their luminosity is smaller the 
larger their core reduction parameter/. This can be seen from 
Fig .121 which shows the presupernova positions of our models in 
the HR diagram. For example, our ‘magnetic’ 16 M© models de¬ 
velop presupernova luminosities around 40 000 L©, whereas the 
16Mo model of Brott et al. (2011) evolves into a presupernova 
red supergiant with a luminosity of 100 000 L©. From Fig.|2]we 
conclude that our ‘magnetic’ models would have to have a mass 
of about 22 M© to produce such a luminous red supergiant pre¬ 
supernova model. 

As expected, models with larger /-parameters develop lower 
core masses (Table [TJ- In fact, the CO cores of our ‘magnetic’ 
16 M© models barely exceed the Chandrasekhar mass. This im¬ 


plies that, within our model assumptions, the limiting initial 
mass for neutron star formation depends on the strength of the 
large scale B field. For the range of core reduction parameters 
explored here, it would obviously not be found far below 16 M©. 
A similar shift may occur for the limiting initial mass for black 
hole formation. 



TeffBK] 



Fig. 2. Evolutionary tracks of stars with M=16, 18 and 20 Mq, com¬ 
puted with different core reduction parameters as indicated, in the HR 
diagram. For the post main sequence phases, small dots are placed on 
the tracks every 105 yr to indicate the speed of evolution of the models 
in the HR diagram. The large circle indicates the presupemova position 
of the models. For comparison, tracks for 15 and 20 Mq stars of non¬ 
rotating models with solar metallicity by Brott et al. (2011) are also 
shown. 

Figure|3] shows the evolutionary tracks of models computed 
with a fixed value of the core reduction parameter/ in the HR 
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diagram, considering the mass range 16Mo to 28 Mq. Whereas 
the 16Mo model evolves into a red supergiant at the beginning 
of core helium burning, the higher mass models all undergo core 
helium burning as blue supergiants. This confirms the trend that 
lower mass models require larger core reduction parameters in 
order to remain on the blue side of the HR diagram. Further¬ 
more, for this fixed value of /, we find no significant trend in the 
effective temperature as a function of mass for our ‘magnetic’ 
core helium-burning blue supergiants. Instead, all the models 
are found in the rather narrow effective temperature range of 
Teff - 17 ... 23 kK. While there seems to be no strong trend with 
mass, this does not imply that all ‘magnetic’ models would pre¬ 
fer this effective temperature range, as indicated, for example by 
our 16 Mq model computed with / = 1.20, which performs core 
helium burning at somewhat lower effective temperatures. 

3.2. Comparison to previous modeis 

While we are not aware of massive star evolutionary models with 
reduced convective cores due to internal large scale magnetic 
fields, models that arise from strong binary interactions can re¬ 
semble our models to some extent. Considering the core helium¬ 
burning stage of evolution, the main feature of our ‘magnetic’ 
models is that the helium core mass is significantly lower than 
that of ordinary stars with the same initial mass. This effect is 
also produced in some close binary models that simulate mass 
transfer, as well as in models for massive star mergers. 

Braun & Langer (1995) simulated binary mass transfer by 
considering single stars that accreted considerable amounts of 
mass on a thermal time scale. They showed that, contrary to 
models using the Schwarzschild criterion for convection, the fi¬ 
nite mixing time obtained when semiconvection is taken into ac¬ 
count may prevent mass gainers from rejuvenating, if the accre¬ 
tion occurs late enough during its core hydrogen-burning phase. 
In such models, the convective core does not adapt to the new, 
higher mass, and thus the helium core mass in the later stages is 
significantly lower than that of ordinary stars. Braun & Langer 
(1995) found that for a given post accretion mass of 20 Mq, 
a blue supergiant core helium-burning evolution was obtained 
when the helium core mass was lower than a certain mass limit. 
This is very analogous to the results we discuss in the previous 
section. 

A similar result is obtained when mass transfer is initiated 
in a massive close binary just after core hydrogen exhaustion 
(Claeys et al. 2011), or in a massive binary merger that occurs 
after the more massive star has exhausted hydrogen in its core 
(Justham et al. 2014). For both situations, core helium-burning 
blue supergiant models are produced. 

4. Comparison with observations 

4.1. Main sequence stars 

For the core hydrogen burning phase of evolution, the prime ob¬ 
servable effect of our assumption that the large scale field leads 
to a lower convective core mass, is that the width of the main 
sequence band in the HR diagram is reduced compared to ordi¬ 
nary models. This effect is clearly visible in Fig.[3j which shows 
the position of the terminal age main sequence (TAMS) of our 
‘magnetic’ models and that of the non-rotating models of Brott 
et al. (2011). 

Figure|3]also shows the positions of the known massive mag¬ 
netic main sequence stars in the HR diagram. These stars are all 
Galactic objects. We find that the observed stars are not equally 



Fig. 3. HR diagram for stars between 16 and 28 Mq with convective 
cores reduced by the same factor, /=1.18. Time steps for every 10^ 
years are shown as small circles, starting at the end of the main se¬ 
quence for clarity. The location of known massive magnetic main se¬ 
quence stars (Briquet et al. 2013; Petit et al. 2013 and references therein; 
Alecian et al. 2014; Fossati et al. 2014; Neiner et al. 2014; Fossati et al. 
2015; Castro et al. 2015) are shown as yellow dots. Additionally, we 
draw the zero age and terminal age main sequences of the non-rotating 
models of Brott et al. (2011) as dashed lines. 

distributed within the main sequence band of the Brott models, 
but that, with some exceptions, they indeed concentrate in the 
area covered by the tracks of our ‘magnetic’ models. We note 
that in a magnitude-limited sample, which roughly holds for the 
observed magnetic stars, one instead expects the cool side of the 
main sequence band to contain more stars than the hot side, since 
cool stars of the same luminosity are visually brighter. This re¬ 
mains true even when considering that for a constant magnetic 
flux, the magnetic held strength is less in more extended stars, 
and the held is thus harder to detect (Fossati et al. 2015). Alter¬ 
natively, the distribution of the magnetic stars in Fig. Ocan also 
be explained by assuming that their B-fleld decays with a time 
scale close to or shorter than their evolutionary time scale (see 
again Fossati et al. 2015). 

While this comparison does not give a final answer, it nev¬ 
ertheless shows that the assumption of magnetic stars having 
smaller convective cores than those of their non-magnetic coun¬ 
terparts leads to a possible interpretation of the data. In this con¬ 
text, it is interesting to mention that the convective core mass in 
massive stars may be measured directly through asteroseismol- 
ogy, as suggested by Mazumdar & Antia (2001) and Mazumdar 
et al. (2006). A comparative investigation of the convective core 
masses of magnetic and non-magnetic stars through asteroseis- 
mology by Briquet et al. (2007, 2012) indeed supports the idea 
that the magnetic stars do develop less massive cores. 

If magnetic main sequence stars are really produced by the 
merging of two main sequence stars in a close binary system 
(Langer 2014), and the large scale held leads to smaller con¬ 
vective cores, then it is difficult to predict to what extent these 
objects will produce blue stragglers in star clusters (Schneider 
et al. 2015). Such stars are more luminous and bluer than the 
bulk of the main sequence stars, meaning that they are located 
beyond the turnoff point in the HR-diagram of the cluster. Even 
if the merged stars successfully rejuvenate, the subsequent sup¬ 
pression of the convective core mass may limit or even cancel 
the blue straggler effect; that is to say the magnetic main se¬ 
quence stars may or may not show up as blue stragglers. On the 
other hand, blue stragglers may also form from slower accretion 
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events, which are less likely to give rise to large scale magnetic 
fields (Langer 2014). Thus, not all blue stragglers are expected 
to be magnetic stars. 

The observational evidence for magnetic fields in blue strag¬ 
glers is scarce, only about a dozen have been analysed so far. 
For roughly half of them magnetic fields were unambiguously 
detected (see for example Bagnulo et al. 2006). Other high- 
precision analyses could not confirm universal magnetic fields 
are present in blue stragglers (Fossati et al. 2010). 

4.2. Supergiants 

Figure [3 shows furthermore that our ‘magnetic’ stars populate 
the region on the cool side of the TAMS of the non-magnetic 
models of Brott et al. In fact, while observationally, this region 
appears to be populated with stars (Fitzpatrick & Garmany 1990; 
Hunter et al. 2008; Vink et al. 2010), single star models gener¬ 
ally do not spend any significant amount of time in this region 
(Ekstrom et al. 2012; Chieffi and Limongi 2013), such that the 
evolutionary history of the blue supergiants is in fact unclear at 
present (McEvoy et al. 2015). 

Within the scope of the VLT-ELAMES Survey of Massive 
Stars, Evans et al. (2008), Hunter et al. (2008) and Vink et al. 
(2010) all find a steep drop in rotation rates at logg ^ 3.2, 
as well as an increase in the surface nitrogen abundance. The 
quoted gravity defines the TAMS of the Brott et al. (2011) mod¬ 
els around 15 Mq. Both features have been recently confirmed 
within the VLT-ELAMES Tarantula Survey (Evans et al. 2011) 
by McEvoy et al. (2015). Our ‘magnetic’ models do provide an 
explanation for the observed population with logg < 3.2. 

If the field in the magnetic main sequence stars does not de¬ 
cay, our models would then predict that the blue supergiants with 
logg < 3.2 should exhibit a large scale magnetic field. Because 
their radii are larger by a factor of 3... 10 and the main sequence 
stars show field strength of several hundred GauB (Petit et al. 
2013; Eossati et al. 2015), flux conservation would lead to field 
strengths of the order of tens of GauB or less (for these stars see 
Fig.0. These fields might thus be difficult to detect (Shultz et 
al. 2014). 

As mentioned in Sect. l3.2l close binary evolution is also able 
to produce long-lived core helium-burning blue supergiants just 
adjacent to the TAMS. It would thus be interesting to investi¬ 
gate the binary fraction of the observed blue supergiants (Langer 
2012). If they were produced by mass transfer or through a bi¬ 
nary merger, detectable companion stars would hardly be ex¬ 
pected (de Mink et al. 2014). The same would be true in the 
magnetic scenario if the large scale fields were really produced 
by strong binary interaction (cf.. Sect. 1). 

Einally, as shown by our models (Eig. 0, some magnetic 
stars with smaller cores may also evolve into red supergiants. 
We see two ways these might be observationally identified. One 
would be to look for the magnetic field itself. While flux con¬ 
servation would lead to inhibitively small average surface field 
strengths, convection in the deep convective envelope might re¬ 
organise the field near the surface such that locally a much higher 
field strength could be possible. On the other hand, the deep con¬ 
vective envelopes of red supergiants may initiate a dynamo pro¬ 
cess that could lead to fields in all red supergiants (Dorch 2004). 
Indeed, Grunhut et al. (2010) find B fields in a third of the 30 
massive cool supergiants that they investigated. 

The other characteristic of our ‘magnetic’ red supergiant 
models is their low luminosity-to-mass ratio. This could show 
up in the photometry of young star clusters showing red super¬ 
giants which are dimmer than expected from the main sequence 


turn-off. Since the luminosity-to-mass ratio is thought to be pro¬ 
portional to the pulsation period in red supergiants (Gough et al. 
1965; Heger et al. 1997), objects with low luminosity-to-mass 
ratios might also be identified through their pulsation properties. 

4.3. Supernovae 

When SN 1987A exploded in the Large Magellanic Cloud, it was 
a big surprise that the progenitor star turned out to be a blue su¬ 
pergiant (Walborn et al. 1989). While single-star models were 
produced soon thereafter that did predict blue supergiant progen¬ 
itors (e.g. Woosley 1988; Saio et al. 1988; Langer et al. 1989), 
single-star models that include the OPAL opacities (Iglesias & 
Rogers 1996) do not appear to produce blue supergiant presu¬ 
pernova stars at the metallicity of the EMC. 

As shown in Sect.[3 many of our ‘magnetic’ models produce 
blue supergiant supernova progenitors. On the other hand, as dis¬ 
cussed in Sect. 13.21 binary mass transfer and post main sequence 
merger may also produce blue presupernova stars. In this con¬ 
text it is interesting to note that other SN 1987A-like supernovae 
have been observed (e.g., Kleiser et al. 2011, Pastorello et al. 
2012), and in fact Smartt et al. (2009) conclude that up to 3% of 
all supernovae could be of this type, meaning they have a blue 
progenitor star. Convective core suppression by large scale mag¬ 
netic fields may help to explain such a large number of events. 

4.4. Magnetars 

Assuming flux freezing, the B field of the massive magnetic 
main sequence stars compressed to neutron star dimension 
would lead to the strength of magnetar fields (Eerrario & Wick- 
ramasinghe 2006). Indeed, if a massive main sequence stars had 
a B field of only ~ lO^G in their core, the resulting B field in 
the neutron star would be of the order of 10^‘^G, which is two 
orders of magnitude greater than typical neutron star magnetic 
field strengths (Langer 2014). In contrast to the scenario by Dun¬ 
can & Thompson (1992), where the magnetar field forms from 
an extremely rapidly rotating collapsing iron core, magnetars as 
successors of magnetic main sequence stars would form slowly 
rotating neutron stars. 

In the view of the idea that the magnetic fields in massive 
main sequence stars might be produced by strong binary interac¬ 
tion, Clark et al. (2014) propose that the blue hypergiant WDl-5 
was the companion star of the magnetar J1647-45 in the young 
Galactic cluster Westerlund 1. To explain the unusual properties 
of WDl-5, in particular its high helium and carbon surface abun¬ 
dances, they propose that reverse mass transfer from the magne¬ 
tar progenitor to WDl-5 must have occurred. This implies that 
WDl-5 was the initially more massive star and the original mass 
donor in the binary system, and that a strong accretion event 
could have formed a stable large scale field in the mass gainer. 
Clearly, a suppression of the convective core mass in the magne¬ 
tar progenitor would have helped to obtain a supernova reversal, 
i.e. to have the initially less massive star explode first and pro¬ 
duce the magnetar in this event. 

5. Discussion and conclusions 

We pursue the assumption that internal large scale magnetic 
fields in massive stars are capable of reducing the mass of their 
convective cores during central hydrogen burning. While we are 
unable to quantitatively predict this effect, we take a parametric 
approach by assuming that effectively larger temperature gra- 
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clients are needed to produce convection compared to the limit 
set by the Ledoux criterion (Gough & Tayler 1986). Indeed, the 
ability of magnetic fields to restrict convection is documented 
well for low-mass stars. The finding of Jackson et al. (2009) 
that the radii for low-mass M dwarfs in the young open clus¬ 
ter NGC 2516 are up to 50% larger than expected is reconciled 
well by MacDonald & Mullan (2013), who apply the Gough & 
Tayler criterion. Feiden & Chaboyer (2013, 2014), who model 
magnetic low-mass stars based on the technique developed by 
Lydon & Sofia (1985), confirm very similar behaviour by com¬ 
paring them to detached eclipsing binaries. 

We followed the evolution of stars in the mass range 16-28 
Mq from the zero age main sequence roughly until neon ignition. 
Our ansatz allowed us to smoothly vary the extent of the core 
mass reduction and study its effects on the star’s main sequence 
and post-main sequence evolution. 

We showed that the stellar models obtained in this way 
predict various observational consequences. Their lower core 
masses compared to ordinary models lead to an earlier termi¬ 
nation of the main sequence evolution, and to blue supergiant 
core helium-burning stars adjacent to the terminal age main se¬ 
quence of ordinary stars. Both effects appear to be supported by 
observations, because the magnetic massive main sequence stars 
appear to prefer the hot side of the main sequence band, and 
the observed hot blue supergiants are not predicted by single¬ 
star evolution otherwise. Our models may also help to explain 
SN1987A-like supernovae, perhaps more abundantly than binary 
evolution. 

On the other hand, binary evolution per se could produce the 
phenomenon of large scale stable B fields in a fraction of the 
massive main sequence stars, such as through stellar mergers. 
Multidimensional MHD calculations are required to investigate 
this question, and binary evolution models could then assess the 
effects of B fields on the further evolution in more detail. 
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